ABSTRACT: The objective of this research was to estimate genetic parameters and genetic trends (GT) for 305-day milk yield (MY305) and 305-day fat yield (FY305) of purebred Dairy Gir animals of the National Dairy Gir Breeding Program. The restricted maximum likelihood method was used in an animal model. GT were obtained via linear regression and divided into two periods (1935-1992 and 1993-2013 for PL305; 1935-1992 and 1993-2010 for MY305). The estimated heritabilities were 0.23 (MY305) and 0.10 (FY305). The GT (kg/year) values for MY305 in the 2nd period for measured females (25.49), females (26.11), and males (35.13) were higher than those found in the 1st period (2.52; 2.06, and 1.00, respectively). The heritability estimated for MY305 confirmed the possibility of genetic improvement by selection and indicated a lower additive genetic effect on FY305 of purebred animals. The genetic progress for MY305 in all purebred population is denoted by the more expressive gains found from 1990’s, when the first bull catalogs were published.

RESUMO: Objetivou-se estimar os parâmetros genéticos e tendências genéticas (GT) para produção de leite (MY305) e produção de gordura (FY305), ambas em 305 dias, de animais puros Gir Leiteiro, integrantes do Programa Nacional de Melhoramento do Gir Leiteiro. Foi utilizada a metodologia da máxima verossimilhança restrita em modelo animal. As GT foram obtidas via regressão linear e divididas em dois períodos (1935-1992 e 1993-2013 para PL305; 1935-1992 e 1993-2010 para MY305). As herdabilidades foram de 0,23 (MY305) e 0,10 (FY305). Para PL305, as GT (kg/ano) do 2º período para fêmeas mensuradas (25,49), fêmeas (26,11) e, machos (35,13) foram claramente superiores às do 1º período (2,52; 2,06 e 1,00; respectivamente). A estimativa de herdabilidade para MY305 reafirma ser possível melhoramento genético por meio de seleção, enquanto para FY305 sugere uma menor influência genética aditiva em animais puros. O progresso genético para MY305 em toda a população pura está evidenciado pelos ganhos mais expressivos, observados a partir da década de 90, quando foram divulgados os primeiros sumários de touros.

ABSTRACT The objective of this work is to estimate genetic parameters and breeding values to improve embryo and oocyte production, using repeatability and random regression models (RRM) for Gir dairy cattle. We used 11,398 records of ovum pick-up from 1,747 dairy Gir donors and evaluated sixteen different models: the traditional repeatability model and fifteen RRM, each of which considered a different combination of Legendre polynomial regressors to describe the additive genetic and permanent environment effects. The 4G1P model (four regressors for the genetic effect and one regressor for the permanent environment effect) is the most suitable model to analyze the number of viable and total oocytes, while the 3G1P is the best model to analyze the number of cleaved and viable embryos, according to the values of the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). The heritability estimated with the RRM was higher than that estimated with the repeatability model. The high repeatability reported for oocyte and embryo count traits indicates that donors, which had high oocyte and embryo counts in the first ovum pick-up, should maintain this result in the next ovum pick-up. Genetic correlations between adjacent ages were high and positive, while genetic correlations between extreme ages were weak. We observed a reranking of the top sires and females (heifers and cows) over the period evaluated. The reliability of the estimated breeding values by RRM showed changes across age, and the expected genetic gains by RRM are larger. This shows that RRM is most suitable alternative for the evaluation and selection of oocyte and embryo count traits.

O objetivo neste estudo foi avaliar a estrutura genética da população de bovinos da raça Girolando no Brasil. Analisou-se o arquivo de pedigree de 26.969 animais, composto de 3.031 machos e 23.938 fêmeas. O nível de conteúdo de informação do pedigree na geração atual foi 61%, mostrando ser de qualidade moderada. O coeficiente de endogamia médio e o coeficiente de relação médio da população Girolando foram 0,11 e 0,13%, respectivamente. O tamanho efetivo da população, considerando a geração completa traçada, foi 188, acima do nível crítico. Do total de 9.457 ancestrais que contribuíram para a população de referência, 457 explicaram 50% da variabilidade genética da população. O número efetivo de fundadores foi 551 e o de ancestrais 393. O intervalo médio de geração foi de 5,26 anos, sendo ligeiramente maior nas trilhas gaméticas mãe-filho e pai-filha. A partir dos coeficientes estimados, pode-se concluir que a endogamia nos rebanhos da raça Girolando foi de pequena magnitude e que as práticas de acasalamento foram adequadas durante o período avaliado. No entanto, é importante continuar com o monitoramento desses coeficientes a fim de prevenir perda de variabilidade genética

The aim of this study was to evaluate the population structure of Girolando cattle in Brazil. The pedigree file contained 26,969 individuals, from which 3,031 were males and 23,938 were females. The average level of completeness of the pedigree in the current generation was of reasonable quality (61%). Inbreeding and average relatedness coefficients were low: 0.11 and 0.13%, respectively. Estimates of effective population size considering the full generations traced was 188, which is above the critical level range. The number of ancestors that contributed to the reference population was 9,457 animals, from which 457 explained 50% of the genetic variability of the population. The effective number of founders and the effective number of ancestors in this population were, respectively, 551 and 393. The average generation interval was 5.26 years, slightly higher in genetic pathways dam-son and sire-daughter. The inbreeding in the Girolando breed was of small magnitude, indicating that the current practices of mating were adequate during the study period. However, it is important to continue monitoring these coefficients in order to prevent loss of genetic variability

Heritability and genetic correlations for honey (HP) and propolis production (PP), hygienic behavior (HB), syrupcollection rate (SCR) and percentage of mites on adult bees (PMAB) of a population of Africanized honeybees were estimated. Data from 110 queen bees over three generations were evaluated. Single and multi-trait models were analyzed by Bayesian Inference using MTGSAM. The localization of the hive was significant for SCR and HB and highly significant for PP. Season-year was highly significant only for SCR. The number of frames with bees was significant for HP and PP, including SCR. The heritability estimates were 0.16 for HP, 0.23 for SCR, 0.52 for HB, 0.66 for PP, and 0.13 for PMAB. The genetic correlations were positive among productive traits (PP, HP and SCR) and negative between productive traits and HB, except between PP and HB. Genetic correlations between PMAB and other traits, in general, were negative, except with PP. The study permitted to identify honeybees for improved propolis and honey production. Hygienic behavior may be improved as a consequence of selecting for improved propolis production. The rate of syrup consumption and propolis production may be included in a selection index to enhance honeybee traits.

Objetivou-se obter os índices de diversidade genética para o SNP (single nucleotide polymorphism) do íntron 4 do gene da osteopontina (OPN) para 434 animais (87 touros e 347 vacas) participantes do Teste de Progênie da raça girolando no Brasil. Para a amplificação, foram utilizados primers descritos para a raça holandesa, e a diferenciação dos alelos C/T desse SNP foi obtida por meio da técnica de PCR-RFLP. As frequências genotípicas TT (52,53%), CT (38,71%) e CC (8,76%) e as frequências alélicas de T (71,9%) e C (28,1%) indicam que a população encontra-se em Equilíbrio de Hardy-Weinberg (EHW). Apesar de o loco do gene OPN estar em EHW, a frequência superior do alelo T do SNP nesses animais pode sugerir uma tendência de fixação do alelo T na raça. Não foi observada diferenciação entre o grupo de touros e vacas (F ST = -0,018), corroborando a estimativa de equilíbrio da população. Considerando os valores estimados pelo F IS (0,043), é possível que ocorram altos números de indivíduos homozigotos para o alelo T observados na população, em virtude da provável herança desse alelo vindo da raça zebuína, e não a endogamia. Assim, para melhor caracterização do polimorfismo do gene OPN, devem ser realizadas avaliações em maior número de animais, uma vez que só foram avaliados animais participantes do teste de progênie.

The objective was to obtain the indices of genetic diversity for the SNP (single nucleotide polymorphism) of the 4 intron osteopontin gene (OPN) for 434 animals (87 bulls and 347 cows) participants in the Teste de Progênie da raça girolando (Girolando Progeny Test) in Brazil. For amplification, primers used were described for the Holstein breed, and differentiation of alleles C/T SNP that was obtained by PCR-RFLP. Genotype frequencies of TT (52.53%), CT (38.71%) and CC (8.76%) and allele frequencies of T (71.9%) and C (28.1%) indicate that the population is in Hardy-Weinberg principle (HWP). Although the OPN gene locus is in HWP, the higher frequency of allele T of SNP in these animals may suggest a setting-trend of allele T in the race. No difference was observed between the group of bulls and cows (F ST = -0.018), supporting the estimate of population balance. Considering the values estimated by the F IS (0.043), it is likely that high numbers of individuals homozygous for the T allele observed in the population occur because of possible inheritance of this allele coming from the zebu breed, rather than inbreeding. Thus, to better characterize the OPN gene polymorphism, assessments in a larger number of animals must be performed, since only animals that participated in the Progeny Test were assessed.

Os registros de produção de leite no dia do controle das três primeiras lactações de 25,5 mil vacas da raça Holandesa foram utilizados para estimar parâmetros genéticos para produção de leite usando duas alternativas de definição da regressão fixa dos modelos de regressão aleatória (MRA). Os polinômios de Legendre de ordens 4 e 5 foram usados para modelar as regressões das curvas fixas (definidas com base nas médias das produções de leite no dia do controle da população ou de múltiplas sub-populações formadas pelo agrupamento de animais que pariram na mesma idade e estação do ano) e aleatórias (genética aditiva e de ambiente permanente) de lactação. Os critérios de informação de Akaike (AIC) e Bayesiano (BIC) indicaram os modelos que consideraram múltiplas regressões das curvas fixas de lactação como os que melhor se ajustaram aos registros de produção de leite da primeira lactação e os modelos que utilizaram uma única regressão da curva fixa, como os melhores para ajuste das segunda e terceira lactações. As herdabilidades para produção de leite ao longo da lactação não variaram entre modelos, entretanto variaram de 0,22 a 0,34; 0,11 a 0,21 e 0,10 a 0,20, respectivamente, para as três primeiras lactações. Semelhantemente às estimativas de herdabilidade os valores das estimativas de correlações genéticas não variaram entre modelos. O uso de uma ou de múltiplas regressões das curvas fixas de lactação pelos MRA não influencia na estimação de parâmetros genéticos para as produções de leite ao longo das lactações de vacas da raça Holandesa.

Records of test-day milk yields of the first three lactations of 25,500 Holstein cows were used to estimate genetic parameters for milk yield by using two alternatives of definition of fixed regression of the random regression models (RRM). Legendre polynomials of fourth and fifth orders were used to model regression of fixed curve (defined based on averages of the populations or multiple sub-populations formed by grouping animals which calved at the same age and in the same season of the year) or random lactation curves (additive genetic and permanent enviroment). Akaike information criterion (AIC) and Bayesian information criterion (BIC) indicated that the models which used multiple regression of fixed lactation curves of lactation multiple regression model with fixed lactation curves had the best fit for the first lactation test-day milk yields and the models which used a single regression of fixed curve had the best fit for the second and third lactations. Heritability for milk yield during lactation estimates did not vary among models but ranged from 0.22 to 0.34, from 0.11 to 0.21, and from 0.10 to 0.20, respectively, in the first three lactations. Similarly to heridability estimates of genetic correlations did not vary among models. The use of single or multiple fixed regressions for fixed lactation curves by RRM does not influence the estimates of genetic parameters for test-day milk yield across lactations.

Avaliaram-se modelos estatísticos com diferentes pressuposições para definir o que melhor descreva a presença de interação genótipo × ambiente no ganho de peso pós-desmama ajustado (GPD345) de bovinos Hereford, mediante o estudo de normas de reação ao ambiente, obtidas por regressão aleatória, usando uma abordagem bayesiana. Quatro modelos hierárquicos de normas de reação (MHNR) foram empregados por meio do programa INTERGEN. O MHNR K utiliza as soluções de grupos contemporâneos estimadas previamente pelo modelo animal padrão (MA) e as considera como nível ambiental para predizer as normas de reação e o MHNR S, que estima simultaneamente esses dois conjuntos de incógnitas. Para ambos os modelos, foram consideradas duas versões, uma com variância residual homogênea (hm) e outra heterogênea (ht). Pelo critério de informação da deviance e fator de Bayes, o MHRN Shm apresentou melhor ajuste aos dados e, pela deviance baseada na ordenada preditiva condicional, o melhor ajuste foi do MHNR Kht, enquanto, pelos três critérios, o pior ajuste foi obtido pelo modelo animal padrão. As herdabilidades estimadas nos MHNR foram crescentes nos gradientes ambientais em GPD345 de -60; 0 e +60 kg. As correlações genéticas estimadas entre o nível e inclinação das normas de reação foram de alta magnitude, entre 0,97 e 0,99, caracterizando efeito de escala em interação genótipo × ambiente. Os modelos hierárquicos de normas de reação são eficientes para descrever as alterações nos componentes de variância decorrentes do ambiente e para descrever a presença de interação genótipo × ambiente na característica GPD345 em bovinos Hereford.

It was evaluated statistical models with different assumptions to define the one that best describes the presence of genotype × environment interaction on adjusted post-weaning weight gain (PWG345) of Hereford cattle, through the study of reactions norms to the environment, obtained by random regression using a Bayesian approach. Four reaction norms hierarchical models (RNHM) were used through the INTERGEN program. The RNHM K uses the solutions of contemporary groups previously estimated by the standard animal model (AM) and considers them as environmental level for predicting the reaction norms and the RNHM S, which jointly estimate these two sets of unknowns. For both models, two versions were considered, one with a homogeneous (hm) and another with a heterogeneous (ht) residual variance. Based on the deviance information criterion and Bayes factor, RNHMs hm showed the best fit to the data, and by the deviance based on conditional predictive ordinate, the best fit was the RNHM Kht, whereas, by all the three criteria used, the worst fit was obtained by using the standard animal model. Heritabilities estimated on RNHM were increasing in the environmental gradients for PWG345, at -60 kg, 0 and +60 kg. The genetic correlation estimated between the level and slope of reaction norms was high, from 0.97 to 0.99, characterizing a scale effect on genotype × environment interaction. The reaction norms hierarchical models are efficient to describe the changes in variance components due to the environment and to describe the presence of genotype × environment interaction on PWG345 trait of Hereford cattle.

Esta pesquisa foi realizada com o objetivo de estimar parâmetros genéticos para as produções de leite, gordura e proteína do leite de vacas da raça Holandesa usando, respectivamente, 56.508, 35.091 e 8.326 registros leiteiros no dia do controle de 7.015, 4.476 e 1.114 vacas, filhas de 359, 246 e 90 touros. Os efeitos genético aditivo e de ambiente permanente foram estimados por REML. Utilizaram-se modelos de regressão aleatória com polinômios de Legendre de ordem 3 a 6. As variâncias residuais foram consideradas homogêneas ao longo da lactação. As estimativas dos componentes de variância apresentaram tendências semelhantes com o aumento da ordem do polinômio em cada característica. As estimativas de herdabilidade variaram, respectivamente, de 0,14 a 0,31; 0,03 a 0,21; e 0,09 a 0,33 para as produções de leite, de gordura e de proteína. As correlações genéticas entre produções de leite, gordura e proteína do leite variaram de 0,02 a 1,00; 0,34 a 1,00 e 0,42 a 1,00, respectivamente. Os modelos com polinômio de Legendre de maior ordem são os mais adequados para ajuste da produção no dia do controle das três características produtivas.

The objective of this study was to estimate genetic parameters for milk, fat and protein yields of Holstein cows using 56,508; 35,091 and 8,326 test-day milk records from 7,015, 4,476 and 1,114 cows, calves of 359, 246 and 90 bulls, respectively. The additive genetic and permanent environmental effects were estimated using REML. Random regression models with Legendre polynomials from order 3 to 6 were used. Residual variances were considered homogeneous over the lactation period. The estimates of variance components showed similar trends, with an increase of the polynomial order for each trait. The heritability estimates ranged from 0.14 to 0.31; 0.03 to 0.21 and 0.09 to 0.33 for milk, fat and protein yield, respectively. Genetic correlations among milk, fat and protein yields ranged from 0.02 to 1.00; 0.34 to 1.00 and 0.42 to 1.00, respectively. Models with higher order Legendre polynomials are the best suited to adjust test-day data for the three production traits studied.

Os totais de 56.508, 35.091 e 8.326 registros, respectivamente, de produção de leite, de gordura e de proteína no dia do controle foram usados para estimar parâmetros genéticos para seis medidas de persistência na produção de leite, de gordura e de proteína de vacas da raça Holandesa criadas em rebanhos do Estado de Minas Gerais. Os componentes de covariância para os efeitos genético aditivo e de ambiente permanente foram estimados via REML por modelos de regressão aleatória com polinômios de Legendre de ordens 3 a 6. Em geral, os modelos com as mais altas ordens dos polinômios de Legendre apresentaram a melhor qualidade no ajuste desses registros produtivos. As estimativas de herdabilidade obtidas pelos modelos para as persistências nas produções de leite, de gordura e de proteína variaram, respectivamente, de 0,04 a 0,32; 0,00 a 0,23; e 0,00 a 0,27. Os valores das estimativas de correlação genética entre persistência e produções de leite, de gordura e de proteína, em 305 dias, variaram de -0,38 a 0,54; -0,39 a 0,97; e -0,78 a 0,67, respectivamente. A medida de persistência proposta por Jakobsen (PS2) é preferencial para uso em futuras avaliações genéticas para persistências nas produções de leite, de gordura e de proteína de animais da raça Holandesa no Estado de Minas Gerais.

Total numbers of 56,508, 35,091 and 8,326 records of milk, fat, and protein test-day yields, respectively, were used to estimate genetic parameters for six persistency measures on milk, fat and protein productions of Holstein cows reared in Minas Gerais state. Covariance components for additive genetic effects and permanent environmental effects were estimated by REML in random regression models using Legendre polynomials from the third to the sixth order. Overall, models with the highest orders of Legendre polynomials showed the best quality of adjustments of these productive records. Heritability estimates obtained by the models for persistence in milk, fat, and protein yields ranged from 0.04 to 0.32, from 0.00 to 0.23, and from 0.00 to 0.27, respectively. Values of genetic correlation estimates between persistence and total 305-day milk, fat, and protein yields ranged from -0.38 to 0.54, from -0.39 to 0.97, and from -0.78 to 0.67, respectively. Persistence measurement proposed by Jakobsen (PS2) is preferential for using in further genetic evaluations for persistence in milk, fat and protein yields of Holstein cows in Minas Gerais state.

Foram utilizados 21.702 registros de produção de leite no dia do controle de 2.429 vacas primíparas da raça Holandesa, filhas de 233 touros, coletados em 33 rebanhos do Estado do Rio Grande do Sul, entre 1992 e 2003, para estimar parâmetros genéticos, para três medidas de persistência (PS1, PS2 e PS3) e para a produção de leite até 305 dias (P305) de lactação. Os modelos de regressão aleatória ajustados aos controles leiteiros entre o sexto e o 300o dia de lactação incluíram o efeito de rebanho-ano-mês do controle, a idade da vaca ao parto e os parâmetros do polinômio de Legendre de ordem quatro, para modelar a curva média da produção de leite da população e os parâmetros do mesmo polinômio, para modelar os efeitos aleatórios genético-aditivo direto e de ambiente permanente. As estimativas de herdabilidade obtidas foram 0,05, 0,08 e 0,19, respectivamente, para PS1, PS2 e PS3 e 0,25, para P305 sugerindo a possibilidade de ganho genético por meio da seleção para PS3 e para P305. As correlações genéticas entre as três medidas de persistência e P305, variaram de -0,05 a 0,07, indicando serem persistência e produção, características determinadas por grupos de genes diferentes. Assim, consequentemente, a seleção para P305, geralmente praticada, não promove progresso genético para a persistência.

There were used 21,702 test day milk yields from 2,429 first parity Holstein breed cows, daughters of 2,031 dams and 233 sires, distributed over 33 herds in the state of Rio Grande do Sul, from 1992 to 2003. Genetic parameters for three measures of lactation persistency (PS1, PS2 e PS3) and for milk production to 305 days (P305) were evaluated. A random regression model adjusted by fourth order Legendre polynomial was used. The random regression model adjusted to test day between the sixth and the 305th lactation day included the herd-year-season of the test day, the age of the cow at the parturition effects and the order fourth Legendre polinomial parameters, for modeling the milk production average curve of the population, and parameters of the same polinomial for modeling the random additive genetic and permanent environmental effects. The estimated heritabilities were 0.05, 0.08 and 0.19, respectively to PS1, PS2 and PS3, and 0.25 to P305, suggesting the possibility of a genetic gain by selection for PS3 and P305. The genetic correlations between persistency measurements and P305 ranged from -0.05 to 0.07, suggesting being, persistency and milk yield, characteristics determined by different gene groups, and that the selection for P305, usually done, do not promote genetic progress for persistency.

Registros da produção de leite de 3.531 controles de 620 primíparas da raça Jersey foram utilizados para a estimação dos componentes de covariância e dos parâmetros genéticos para a produção de leite por meio dos modelos de repetibilidade e de regressão aleatória. O modelo de repetibilidade (MR) incluiu os efeitos fixos de rebanho-ano-mês do controle, estação do parto e idade da vaca ao parto como covariável, com termos linear e quadrático, e os efeitos aleatórios de animal, de ambiente permanente e residual. O modelo de regressão aleatória incluiu os mesmos efeitos do MR, com o uso da função de Wilmink para a modelagem dos efeitos fixo (estação de parto) e aleatórios (genético de animal e de ambiente permanente). As estimativas de variância genética aditiva, de ambiente permanente, residual e fenotípica, e de herdabilidade e repetibilidade para a produção de leite usando-se o MR foram, respectivamente, 7,08; 1,88; 7,79 e 16,76 kg² e 0,42 e 0,50. Pelo modelo de regressão aleatória, os valores de herdabilidade para a produção de leite decresceram de 0,81 a 0,20 ao longo da lactação. As estimativas de variância genética aditiva mantiveram-se praticamente constantes durante a lactação, com maiores valores no início e final da lactação. As correlações genéticas foram próximas à unidade para as produções de leite da fase intermediária da lactação e diminuíram com o aumento do período entre os controles leiteiros. O modelo de repetibilidade não é a melhor alternativa para a modelagem das produções de leite no dia do controle da primeira lactação de vacas Jersey no Brasil.

Data consisting of 3,531 test day records of 620 first parity Jersey cows were used to estimate genetic parameters for milk yield. Covariance components were estimated by repeatability and random regression animal models. The Wilmink function was used in the random regression model to fit the fixed effect of season of calving and the additive genetic, permanent environment and residual random effects. The repeatability model included the fixed effects of herd-year-month of test, season of calving and age at calving as a covariable with linear and quadratic terms and animal, permanent environmental and residual random effects. The estimates of additive genetic, permanent environmental, residual and phenotypic variances, heritability and repeatability obtained by the repeatability model for milk yield were, respectively: 7.08; 1.88; 7.79 and 16.76 kg², 0.42 and 0.50. The estimates of heritability for milk yield obtained from the random regression model decreased from 0.81 in the beginning of the lactation to 0.20 in the end. Larger estimates of genetic variance were observed at the extremes of the lactation period, but they were practically constant in mid-lactation. The permanent environmental variances increased along lactation. Genetic correlations between milk yields were larger in close test days, particularly in mid-lactation and decreased for the most distant test days. Genetic parameters do not indicate the repeatability model as an alternative to use test day records in order to estimate breeding values for milk yield of Jersey cows in Brazil.

Registros de produção de leite de 68.523 controles leiteiros de 8.536 vacas da raça Holandesa, filhas de 537 reprodutores, distribuídas em 266 rebanhos, com parições nos anos de 1996 a 2001, foram utilizados na comparação de modelos de regressão aleatória, para estimação de componentes de variância. Os modelos de regressão aleatória diferiram entre si pelo grau do polinômio de Legendre utilizado para descrever a trajetória da curva de lactação dos animais. Os modelos incluíram os efeitos rebanho-mês-ano do controle, composição genética dos animais, freqüência de ordenhas diárias, regressão polinomial em cada classe de idade-estação de parto para descrever a parte fixa da lactação e regressão polinomial aleatória relacionadas aos efeitos genético direto e de ambiente permanente. As estimativas de herdabilidade obtidas oscilaram de 0,122 a 0,291. Verificou-se que o modelo de regressão aleatória que utilizou a maior ordem para os polinômios de Legendre descreveu melhor a variação genética da produção de leite, de acordo com o critério de Akaike.

Data comprising 68,523 test day milk yield of 8,536 cows of the Holstein breed, daughters of 537 sires, distributed in 266 herds, calving from 1996 to 2001, were used to compare random regression models, for estimating variance. Test day records (TD) were analyzed by different random regression models regarding the function used to describe the trajectory of the lactation curve of the animals. Legendre orthogonal polynomials function of second, third and fourth order were used. The random regression models included the effects of herd-month-year of the control, genetic group of the animals; the frequency of the daily milk; regression coefficients for each class of age-season (in order to describe the fixed part of the lactation curve) and random regression coefficients related to the direct genetic and the permanent environmental effects. The heritability estimates obtained using the random regression models ranged from 0.122 to 0.291. The random regression model which used the fourth order Legendre polynomials was the model which better described the genetic variation of the milk yield, according to AIC test.

Foram estudados os registros da idade ao primeiro parto (IPP), intervalo de partos (IDP) e período de serviço (PSERV) provenientes das cinco primeiras lactações de 1.406 vacas do ecótipo Mantiqueira, filhas de 113 reprodutores, com partos entre os anos de 1977 e 1997, do Arquivo de Escrituração Zootécnica do Pólo Regional do Vale do Paraíba da Agência Paulista de Tecnologia dos Agronegócios da Secretaria de Agricultura e Abastecimento do Estado de São Paulo (NPZGJRA/IZ/APTA/SAA-SP), no município de Pindamonhangaba. Os modelos usados para obtenção das médias dos quadrados mínimos incluíram, além dos efeitos aleatórios de reprodutor e do erro, os efeitos fixos de ano e estação de nascimento (IPP) ou ano e estação de parto (IDP e PSERV) e idade ao parto. Para a obtenção dos parâmetros genéticos pela metodologia REML, em análises de uni e bicaracterísticas, foram utilizados modelos que ajustavam os efeitos fixos de ano-estação e idade da vaca ao parto, além dos efeitos aleatórios de animal e do erro. As médias ajustadas ± erros-padrão foram: 45,60 ± 8,14 meses para a IPP; 460,56 ± 122,05 dias para o IDP; e 181,44 ± 110,31 dias para o PSERV. As herdabilidades foram: 0,13; 0,01; e 0,01, respectivamente, para IPP, IDP e PSERV.

The reproductive performance of Mantiqueira cows were analyzed using data of records of age of first calving (IPP), calving interval (IDP) and days open (PSERV). These data were the first five lactations of 1406 cows of Mantiqueira ecotype consisting of daughters of 113 sires, with calving from 1952 to 1997, from the Paraiba Valley Regional Center of Research in São Paulo State (APTA/SAA-SP) breeding program. Models to obtain least square means included fixed effects for year and season of birth (for IPP) or calving (for IDP and PSERV) and age, in addition to the random effect of sire and error. The models to obtain genetic parameters by REML methodology included fixed effects for year-season and age, in addition to the random effects of animal and error. The least square means ± standard error were 45.60 ± 8.14 months to the IPP, 460.56 ± 122.05 days to the IDP and 181.44 ± 110.31 days to the PSERV. Estimated values of heritabilities were: 0.13 (IPP), 0.01 (IDP), and 0.01 (PSERV).